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We investigate trilepton final states to probe top anomalous couplings at the Large Hadron Col- 
lider. We focus on events originating from the associated production of a single top quark with a 
Z-boson, a channel sensitive to several flavor-changing neutral interactions of top and up/charm 
quarks. In particular, we explore a way to access simultaneously their anomalous couplings to Z- 
bosons and gluons and derive the discovery potential of trilepton final states to such interactions 
with 20 fb _1 of 8 TeV collisions. We show that effective coupling strengths of 0(0.1 - 1) TeV" 1 
can be reached. Equivalently, branching fractions of top quarks into lighter quarks and gluons or 
Z-bosons can be constrained to be below 0(0.1 — 1)%. 
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INTRODUCTION 



Since its discovery, the top quark, given its large 
mass, is generally considered as a sensitive probe to 
new physics. In particular, its flavor-changing neutral 
couplings to gluons and Z-bosons gqt and Zqt (with 
q = u,c) vanish at tree-level as consequences of the unbro- 
ken QCD gauge symmetry and the Glashow-Iliopoulos- 
Maiani mechanism. This also guarantees that in the 
Standard Model, these interactions stay suppressed at 
higher orders. Deviations from these predictions have 
therefore been widely searched for at hadron colliders in 
flavor-changing neutral top decays |l|— 15| and single top 
processes p-lpj. Expressed in terms of branching ratios, 
current constraints are given at the 95% confidence level 
by BTl{t^ug)<b.l ■ 1CT 3 %, BK(t ->• eg) < 2.1 • 1(T 2 % 
and B7Z(t — > qZ) < 0.21%, assuming one single non- 
vanishing anomalous coupling. 

Existing phenomenological analyses J10NL5| mostly fo- 
cus on direct top quark production and, in a smaller 
extent, on rarer processes and indirect probes. In this 
work, we concentrate on the associated production of a 
top quark with a Z-boson. Although monotop signa- 
tures related to invisible Z-boson decays can be interest- 
ing for new physics searches [16| , we focus on final states 
where both heavy particles decay into charged leptons. 
On the one hand, this yields a trilepton topology infer- 
ring a Standard Model background under good control. 
On the other hand, this allows us to design an analysis 
probing at the same time gqt and Zqt interactions, in con- 
trast to standard investigations that are only sensitive to 
one type of couplings. Since the production rate of a top 
quark in association with a Z-boson is highly suppressed 
in the Standard Model, we aim to interpret possible ex- 



cesses in trilepton events as hints of large anomalous gqt 
and Zqt couplings, as predicted by many new physics 
theories [13]. 

In order to facilitate the analysis of new physics effects 
in top anomalous couplings, the relevant interactions are 
usually described in terms of a minimal set of effective 
operators independent of the underlying theory [18| . This 
bottom-up approach is in particular well motivated when 
new particles are heavy so that they can be integrated 
out. In this context, limits on new physics can be derived 
in a modcl-indcpcndcnt way. To this aim, we analyze, by 
means of Monte Carlo simulations, trilepton events to be 
produced at the Large Hadron Collider (LHC). 

This work extends a pioneering study focusing on 
14 TeV runs of the LHC with respective luminosity of 
10 ftr 1 and 100 fb" 1 (nj- First, we update the older 
results according to current LHC settings, i.e., 20 fb _1 
of proton-proton collisions with a center-of-mass energy 
of 8 TeV. Next, we make use of a more accurate de- 
scription of the Standard Model background relying on 
state-of-the-art Monte Carlo event generation including 
multiparton matrix-clement merging and on a more ad- 
vanced simulation of the detector response. We then de- 
sign a better adapted search strategy that we believe to 
be worthwhile to be tested by the ATLAS and CMS ex- 
periments. 



EVENT SIMULATION 

Single top production in association with a Z-boson 
can be driven, at tree-level, by two main mechanisms re- 
lated to the strong and weak sector, respectively. This is 
illustrated by the two representative Feynman diagrams 
of Figure [1] where flavor-changing interactions of top and 




FIG. 1: Representative Feynman diagrams leading to the pro- 
duction of a single top quark in association with a Z-boson. 
Flavor violation (red dots) occurs either in the strong (left) 
or weak (right) sector. 



light 1 quark fields t and q involve either a gluon (G") or a 
Z-boson (Zfj). Both cases can be described by means of 
effective operators to be supplemented to the Standard 



u 



Model Lagrangian 



19] . Denoting the fundamen- 
tal representation matrices of SU(3) by T a , the strong 
and weak coupling constants by g s and g, the cosine of 
the weak mixing angle by cw and the Z-boson and gluon 
field strength tensors by Z M „ and G" , this effective La- 
grangian reads 
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In this expression, Pj, and P R are chirality projectors 
acting on spin space and a^ v = |[7 M ,7 l/ ]. While other 
operators are in principle possible, they can always be re- 
cxprcssed in terms of those included in Eq. (Q} so that we 
only consider the minimal set of independent Lagrangian 
terms above. The magnitude of new physics effects that 
are assumed to appear at an energy scale A are mod- 
eled through real dimensionless parameters n gqt , n zq t and 
Czqt together with complex chiral parameters f q ' , f q ' R 
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\f q Y = 



and f q > R normalized to \f q \ 2 + \f, 

l/^| 2 + l/f| 2 = l. 

In this work, we perform a phenomenological analy- 
sis based on Monte Carlo simulations of collisions pro- 
duced at the LHC, running at a center-of-mass energy of 
■s/s =8 TeV and for an integrated luminosity of 20 fb _1 . 
For both sig nal and background, wc use the Mad- 
Graph 5 [20j package for the generation of hard scatter- 
ing matrix elements including up to two additional jets 
and convolved with the leading order set of the CTEQ6 
parton density fit [21] . Parton- level events produced 
in this way are then matched to parton showering and 
hadronization by means of the program Pythia 6 12211 
and merged according to the fc^-MLM scheme 23J, |24j . 



Detector effects arc subsequently accounted for using a 
modified version of Delphes 2.0 [25]. The latter includes 
a modeling of the performances of the CMS detector as 
described in Ref. [261 ] and a more recent description of 
the 6-tagging efficiency and mistagging rates. For the 
latter, we base our implementation on the Tchel al- 
gorithm [27|, [28] , which leads, e.g., to a ~ 75% tagging 
efficiency and a 10 — 15% mistagging rate for jets with 
a transverse momentum of 50 — 80 GeV. After employ- 
ing the Fast Jet package [29j for jet reconstruction with 
an anti-/c t algorithm of radius parameter set to R = 0.4, 
simulated events are analyzed within the MadAnaly- 
SIS 5 framework J30j . 

In order to allow for signal simulation, we have im- 
plemented the effective Lagrangian of Eq. ([TJ within the 
FeynRules package 3ll435] and subsequently exported 
the model to a UFO module [36] that has been linked 
to MadGraph 5. We then generate events describing 
the production of a top (anti) quark decaying lcptonically 
(together with a 6- jet and missing energy), in association 
with a pair of same flavor leptons with opposite elec- 
tric charges 2 whose the invariant mass is greater than 
12 GeV. In contrast to the diagrams of Figure [TJ virtual 
photon and off-shell Z-boson effects are included in the 
simulation. 

In our study, we investigate simplified scenarios with 



fq = # = fq = ° and fq = /* = /* = L Leading order 



inclusive cross sections a t z{x) read, assuming a single 
non-vanishing anomalous coupling x at a time, 



(Jtz{K gut /K) = 8.678 



crtz(«gct/A) = 3.255 
<?tz(Kzut/A) = 5.769 

o-«j(«W A ) = 2 - 729 
<rtz(<*«t/A) = 2.727 
<rtz(Czct/A)= 1-533 
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Those predictions agree with the results of Ref. |37j and 
motivate us to neglect possible effects of the ( zq t pa- 
rameters that lead to significantly lower cross sections, 
not reachable within 2012 LHC data, for anomalous cou- 
plings of O(0.1 — 1) TeV -1 . In addition, we also consider 
that strong (n gq t) and weak (n zq t) contributions can be 
treated independently and we include next-to-lcading or- 
der (NLO) AT-factors fixed to 1.3 Q. 

We now turn to the simulation of the Standard Model 
background. First, we do not consider multijet events 
since their correct treatment requires data-driven meth- 
ods. We instead rely on existing experimental analy- 



1 The naming light quark (jet) also refers to charmed objects. 



2 By the terminology leptons, we equivalently denote electrons, 
muons and taus decaying lcptonically 



ses that have shown that they contribute negligibly after 
a selection based on the trilepton topology and missing 
energy requirements [38], |39(. Events originating from 
the production of a single gauge boson, decaying lepton- 
ically or invisibly, together with jets have been normal- 
ized to the next-to-next-to- leading order (NNLO), using 
total rates of 35678 pb and 10319 pb for W- and Z- 
boson production calculated by means of the Fewz pro- 
gram 40j, |41| with the CT10 parton densities 42j. Inclu- 



sive top-antitop events have been normalized to 255.8 pb, 
as predicted by the Hathor package [43] after convolv- 
ing all NLO diagrams and genuine NNLO contributions 
with the CT10 parton densities. Single top event genera- 
tion has been split into the production of three inclusive 
samples, normalized, at an approximate NNLO accuracy, 
to 87.2 pb, 22.2 pb and 5.5 pb for the t, tW and s chan- 
nels [44|. Diboson event weights have been rescaled ac- 
cording to NLO predictions of 30.2 pb, 11.8 pb and 4.5 pb 
for the WW, WZ and ZZ modes, respectively, as com- 
puted by means of the Mcfm program [45j-|47[ , employing 
again CT10 densities and after neglecting full hadronic 
decay modes. Finally, ttV (with V — W,Z) events have 
been normalized to NLO as predicted by Mcfm while 
other simulated rare Standard Model processes rely on 
MadGraph 5 leading-order results. Inclusive cross sec- 
tions for the ttW, ttZ, tZj, ttWW and tttt channels have 
been found to be 0.25 pb, 0.21 pb, 46 fb, 13 fb and 0.7 fb, 
respectively. 



ANOMALOUS TOP COUPLINGS AT THE LHC 

Events are preselected by requiring exactly three iso- 
lated charged leptons (electrons or muons) with a trans- 
verse momentum px > 20 GeV and a pseudorapidity 
1 77 1 < 2.5, using the CMS coordinate system presented, 
e.g., in Ref. 26|. Lepton isolation is enforced by opti- 
mizing, independently for electrons and muons, an isola- 
tion variable I IC \ so that the analysis sensitivity to top 
anomalous couplings is maximized, the sensitivity be- 



ing defined as the significance S/vS + B where S and 
B are respectively the number of signal and background 
events after all selections. The quantity / ro i is derived 
from the amount of transverse energy, evaluated rel- 
atively to the lepton px, present in a cone of radius 
R= \/Atp 2 + A77 2 =0.3 centered on the lepton, ip being 
the azimuthal angle with respect to the beam direction. 
The selected electrons and muons are asked to satisfy the 
constraint I IC \ < 0.28 and I IC \ < 0.06, respectively. 

Two leptons, labeled by l\ and £2, are tagged as origi- 
nating from a Z-boson decay by selecting the pair of same 
flavor leptons with opposite electric charges whose the in- 
variant mass rri£ t i 2 is the closest to the Z-boson mass mz- 
We further reject events for which \mi 1 i 2 —mz\> 15 GeV. 

The third lepton £3 is assumed to originate from a lep- 
tonically decaying top quark, being thus accompanied by 



TABLE I: Number of expected events for 20 fb" 1 of LHC 
collisions at ,/s — 8 TeV, together with the associated sta- 
tistical uncertainties, after each of the selections described in 
the text. Two representative signal scenarios have been con- 
sidered, with non-vanishing n gu t (second column) and k zu j 
(third column) parameters, respectively, after enforcing a sig- 
nal normalization of 10 fb. The sum over all background 
contributions is also indicated (last column). 



Selection 


Kgut 7^ 


Kzut -£ 


Background 


Trilepton topology 


41.3L0.3 


36.6L0.3 


3648.3L143.1 


m ilia € [75, 105] GeV 


38.1L0.3 


34.0L0.3 


3271.8L136.0 


fir > 30 GeV 


31.8L0.2 


26.9L0.3 


1362.9L31.9 


rn^ > 10 GeV 


29.6L0.2 


25.4L0.2 


1256.9L18.8 


At least one jet 


25.2L0.2 


22.4L0.2 


549.0L18.5 


One single fa-tagged jet 


17.4L0.2 


14.7L0.2 


102.8L2.3 


m (reco) < 25Q ^y 


16.4±0.2 


13.8±0.2 


55.2L1.9 



a neutrino yielding missing transverse energy fix and by 
a jet issued from the hadronization of a 6-quark. We 
therefore select events with fix > 30 GeV and require the 
reconstructed VF-boson transverse mass m^ (as defined, 
e.g., in Ref. (5j) to fulfill m^ > 10 GeV. In addition, 
we focus on events containing at least one jet satisfy- 
ing pt > 30 GeV and \q\ < 2.5 and exactly one 6-tagged 
jet. For all jets, the ratio between the hadronic and elec- 
tromagnetic calorimeter deposits is also required to be 
larger than 30%. Finally, the four- momentum of the W- 
boson from which the missing energy and the lepton £3 
are issued is reconstructed, using the VF-boson mass as 
a constraint. The top quark mass mj re can then be 
determined by computing the invariant mass of the W- 
boson and the 6-tagged jet system, and further employed 
for background rejection. The reconstructed top mass 



is asked to fulfill 



(rcco) 



< 250 GeV, this criterion being 



found to maximize the analysis sensitivity to new physics. 
The number of events surviving each of the selection 
criteria is indicated, together with the associated statisti- 
cal uncertainties, in Table fl] for both the Standard Model 
background and two representative signal scenarios for 
which the only non-vanishing anomalous coupling is re- 
spectively Kg U t and K zu t' After all selections, 34.8 ± 0.7 
diboson, 19.0 ± 1.8 it and 1.3 ± 0.1 tZj events remain, 
all other background contributions being negligible. In 
comparison, 16.4 ±0.2 and 13.8 ±0.2 events are expected 
for the two signal scenarios after normalizing both sam- 
ples to a cross section of 10 fb, the difference being due 
to a different acceptance. 

The effect of the top reconstruction is illustrated on 
Figure [2] After applying the selection strategy described 
above but the last step, we first present the m t re spec- 
trum in the Standard Model. The predictions show a 
broad peaky distribution, with a maximum close to the 
top mass. We then superimpose two representative sig- 
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FIG. 2: Distribution of the reconstructed top mass mj r e 
after the selection strategy described in the text, with the 
exception of the final cut, for 20 fb" 1 of LHC collisions 
at yfs =8 TeV. We indicate the dominant diboson (VV), 
top-antitop pair with possibly one or several gauge bosons 
(ti(+X)) and tZj contributions to the Standard Model back- 
ground, as well as two signal spectra for two representative 
new physics scenarios, normalized to 10 fb. 



nal distributions, normalized to 10 fb. They all present 
narrower peaks, the distribution maximum being at the 
top mass. Confronting the shapes of the two signal spec- 
tra also indicates that observing tZ trilepton events does 
not provide any discriminating power among weak and 
strong flavor-changing neutral top interactions. 

The sensitivity of the 2012 LHC run to the four 
top anomalous couplings under consideration, defined as 
S/y/S + B (see above), is presented for several choices of 
the k/A parameters in Figure [3] in the case we assume a 
single non- vanishing coupling at a time. The results are 
subsequently fitted by polynomial functions so that 3<7 
and 5cr discovery ranges are extracted. Hence, anoma- 
lous couplings such as 

Kgut/ A > 0.12 Te^ 1 (0.09 TeV -1 ) , 
Kgct/A > 0.42 TeV" 1 (0.31 TeV" 1 ) , 
K Z ut/A > 0.56 TeV" 1 (0.40 TeV" 1 ) , 
K zct /A > 1.78 TeV" 1 (1.30 TeV" 1 ) , 



(3) 



can be reached at the 5a (3a) level with 20 fb -1 of 8 TeV 
LHC collisions. Using NNLO results for the decay width 
of the top quark in the Standard Model 48j , the limits 



can be translated in terms of constraints on the branching 
fractions of rare top decays, 



BK(t -> gu) < 0.47% (0.25%) , 
BK(t -> gc) < 5.1% (2.8%) , 
BTZ(t -> Zu) < 0.39% (0.20%) , 
BK{t -> Zc) < 3.8% (2.1%) , 

at the 5(7 (3a) levels. The obtained limits on BlZ(t 



(4) 



_ 


I V 1 I 




1 


L 


r • K gm 






/ i 




I ■ K g« 






I 




_ * k 2U1 






1 










t 


- 










_ 


~ \ls = 8 TeV 








-_ 


I L rnt = 20fb- 1 








-_ 


ZZ^Z^ 




rrVi 




r 



10' 1 



1 10 

Kv qt /A[TeV- 1 ] 



FIG. 3: LHC sensitivity to the considered top anomalous 
couplings as a function of the interaction strengths, given as 
S/y/S + B after applying the analysis described in the text. 
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FIG. 4: 3<7 and 5er discovery ranges for top anomalous inter- 
actions in the strong (blue, non- zero K gut /A and K gc t/A) and 
weak (red, non-zero K zut /A and K zc t/A) sector. 



gq/Zc) are found to be not competitive by one or two or- 
ders of magnitude with results extracted from direct top 
production and flavor-changing top decays p, Q . In con- 
trast, the bounds obtained on BlZ(t — > Zu) are promis- 
ing. First, they show that the trilepton analysis which 
we are proposing can be used to confirm (or possibly 
improve) the current constraints extracted from flavor- 
changing top decays in tt events at y/s — 7 TeV 0, [f| 
by means of an independent measurement at a different 
center-of-mass energy an with a different integrated lumi- 
nosity. Next, both channels are statistically independent 
and thus valuable to be combined for getting stronger 
bounds on the Ztu coupling strength. 

We now allow for several non-zero couplings simultane- 
ously, either in the strong sector (non- vanishing K gu t/A 
and Kg C t/A parameters) or in the weak sector (non- 
vanishing K zut /A and K zc t/A parameters). In Figured 



we extract the associated 3<7 and 5a discovery reaches. 
We observe a better sensitivity to flavor-changing inter- 
actions with an up quark than with a charm quark, as 
expected from parton densities, the charm content of the 
proton being suppressed with respect to its up content. 



CONCLUSIONS 

In this letter, we have investigated the discovery po- 
tential of trilcpton final states to anomalous, flavor- 
changing, neutral couplings of the top and light quarks to 
gluon and Z-boson fields. These interactions induce pos- 
sibly significant production rates of such final states via 
the associated production of a leptonically decaying top 
quark with a dilepton pair. We have designed a search 
strategy capable to probe simultaneously gqt and Zqt 
couplings and showed that the 2012 LHC run can reach 
interaction strengths of order 0.1 — 1 TeV -1 . Equiva- 
lcntly, bounds on the associated rare top branching frac- 
tions as low as 0(0.1 — 1)% could be set, possibly improv- 
ing, in particular, the current limits on the t—>Zu decay. 
Our results therefore motivate a further investigation by 
the ATLAS and CMS collaborations in the context of 
real LHC data. 
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